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FIGUR E 

A . 1 . Variation of functions F and F^ (Arrows correspond to q# 

decreasing) . 

(a) Adiabatic. 

(b) Isothermal. 

B. l. Velocity distribution functions at q$ ■ 0 and q$ < 0. 

(a) Waterbag with all particles forward-going. 

(b) Waterbag with forward- and backward-going particles. 

(c) Maxwellian with large dr if t-to-thermal velocity ratio. 

(d) Half-Maxwellian. 

B.2. Variation of N(t) for the velocity distribution functions 
of figure B.l. 
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A Fluid description of pijvsma double-layers 


by 

J. S. Levine and F. W. Crawford 
institute for Plasma Research 
Stanford University, CA 94305, USA 

ABSTRACT 

This paper describes the space-charge double-layer that forms 
between two plasmas with different densities and thermal energies. 

Three progressively more realistic models are treated by fluid theory, 
taking into account four species of particles: electrons and ions 

reflected by the double-layer, and electrons and ions transmitted 
through it. First, the two plasmas are assumed to be cold, and the self- 
consistent potential, electric field and space-charge distributions within 
the double-layer are determined. Second, the effects of thermal velo- 
cities are taken into account for the reflected particles, and the modifi- 
cations to the cold plasma solutions are established. Third, the further 
modifications due to thermal velocities of the transmitted particles are 
examined. The applicability of a one-dimensional fluid description, 
rather than plasma kinetic theory, is discussed. One valuable product 
of this description is the potential difference across the double-layer 
in terms of the parameters of the two plasmas which it separates. A 
useful length parameter is defined characterizing the distance over 
which most of this potential is dropped. Comparisons are then made 
between theoretical predictions, and double-layer potentials and lengths 
deduced from laboratory and space plasma experiments. 
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1. INTRODUCTION 


Double-layers consist of two space-charge layers In close proximity, 
one positively churned and one negatively. They are commonly observed 
In the laboratory as cathode sheaths (Langmuir 19P9» Crawford & Cannara 
I9f ; Prewett & Allen 1970. or 88 constriction sheaths If the vessel con- 
taining a plasma column has a reduction in cross-sectional area (Crawford 
& Freeston 1963; Andersson et a_l. 1*69; Sandahl 1971; Jacobson 6 Eubank 
1973). 

Double-layers can occur that are not controlled by electrodes or by 
the boundary. Laboratory experiments on positive columns In a variety of 
neutral gases, with (Lutsenko et^ «U. 197 ; Torven & Andersson 1978) and 
without (Bablc Torven 1}7^; Armstrong k Torven 197 1 *; Armstrong I97t I 
Levine et al. 1978) magnetic fields, have demonstrated the onset of double- 
layers as the current density is raised. They have also been observed In 
double- (Quon t. V. ing 197 ' ) and triple-plasma (Coakley et^ «U. 1978 ) devices 
as transitions between plasmas of differing characteristics, and in various 
computer simulations of such plasmas (Goertz & Joyce 1975; DeGroot e£ al . 

1 ) 77 . Joyce s* Hubbard 1978 ; Hubbard ca Joyce 1978 ). The experimental evidence 
for double-layers has been reviewed by Torven ( 1978 ). There Is a growing 
body of evidence that double-layers form in the magnetosphere, accelerating 
the high energy electrons that are associated with auroral displays (see, 
e.g. Shawhan et al. 1978 and references therein). Double-layers have also been 
invoked in discussions of solar flares (Alfven sc Carlqvist 1967; Carlqvist 
19* ); Hasan ter Haar 1978), ana to explain how Io modulates Jovian deca- 
metric radiation (Smith w Goertz 1978 ). 

A variety of theoretical models of state-state double-layers have 
been proposed which can be understood qualitatively from figure 1(a): a 

monotonic potential variation occurs over a length that is long on the 
scale of the electron and ion Debye lengths, but small on the scale of 
laboratory or space dimensions. The localization of the potential step, 
and the associated electric field, implies that although neutrality is 
violated within the double-layer, the charge, integrated across the double- 
layer, is zero. The plasmas in Regions 1 and 2 may be characterized by 
different densities, drift velocities and temperatures. The double-layer 
is assumed to be much shorter than the collisional mean free path, so that 


collision effects within it may be Ignored, though they may be Important 
In Regions 1 and ?, The double-layer Is generally treated as a one- 
dimensional structure, thus excluding magnetic field effects, but a two- 
dimensional model in which the double-layer is oblique to the magnetic 
field, has been treated (Swift 197!')* A general review o^ double-layer 
theory has been given by Carlqvist (1978)* 

Several specific models of double-layers that have been studied are 

shown in figure l(b)-(e). In the first of these, Region 1 is replaced by 

an electron-emitting cathode, injecting a flux r of cold electrons, and 

Region 2 is replaced by an ion-emitting anode, injecting a flux f, of 

cold ions, as shown in figure 1(b). The Langmuir (1929) theory of a space- 

charge-limited diode can be applied to the double-layer. In this, Poisson's 

equation is solved consistently with the particle fluxe . The well known 

Langmuir relation, I /I~ = (m /m )^ is derived from the assumption that 

el i e 7 

the electric field vanishes at the emitting surfaces. 

In a model studied by Crawford ta Cannara (19* 1 ), appropriate to a hot 
cathode discharge. Region 1 is replaced by an electron-emitting cathode, 
and Region 2 is taken to be a uniform infinite plasma, as illustrated in 
figure 1(c). Plasma electrons reflected within the sheath are assumed to 
obey Boltzmann's law. Ions entering the double-layer from the plasma are 
taken as monoenergetic. A minimum ion velocity is found which, in the 
limit of a low flux of electrons injected from the cathode, is the Bohm 
( 19 ^ 9 ) sheath condition, v^ > T^/m^ • 

Block (1972) treats the double-layer between two infinite plasmas, 
using fluid theory to include temperature effects, and taking account of 
a population of ions in Region 1, and of electrons in Region 2, that are 
reflected by the potential rtep, as shown in figure 1(d). For the case of 
an infinite potential difference across the double-layer, Block shows that 
the Langmuir condition is obeyed, and that modified Bohm conditions must 
be satisfied for the model to be self-consistent. 

Montgomery * Joyce (l> 9) use kinetic theory to show that a double- 
layer may be constructed as a stationary shock-like solution of the Vlasov 
equation in a system with no current. As illustrated in figure 1(e), they 
assume two streaming populations in Region 1; electrons that are accelerated 
by the potential, and ions of the same density and velocity that are decele- 
rated, but transmitted through the double-layer. A second electron 
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population Is required on the high potential side of the double-layer to 
ensure electrical neutrality. The transmission of the Ions Imposes an 
upper limit on the potential difference equal to the Ion streaming energy. 

Kan ( 19 T ) extends the Montgomery a Joyce (1969) model to show that 
electrostatic shock solutions can exist In u T » plasma under 
conditions of the high latitude plasma sheet. Populations of reflected 
ions In Hegion 1 and transmitted Ions from Region 2 are Included, and the 
assumption of zero current Is eliminated. Magnetic field effects are 
treated In this one-dlmenslonal theory by using loss-cone distribution 
functions. Necessary conditions for existence of the shock are derived 
In terms of the ratios of drift to thermal velocity for Ions and electrons, 
for an assumed potential difference. 

Andrews Allen (1971) describe double-layers In terms of averages 
over unspecified distribution functions, and use Maxwellian distributions 
for the reflected particles, as In figure 1(d). Numerical results for 
transmitted particles with delta-function distributions are presented, 
which satisfy Bohn conditions, but give values of I ^/T. that are less 
than the Langmuir condition. 

Hasan it ter Haar (I 978 ) analyze the double-layer for delta-function 
and power law distributions for the transmitted particles In figure 1(d), 
using waterbag distributions for the reflected particles. Conditions 
analogous to the Bohm conditions are derived, and a modified Langmuir con- 
dition is derived for delta-function distributions of transmitted particles. 

Using a different approach, Knorr it Coertz (197*0 assume forms for 
the potential and for the velocity distributions of three of the four particle 
species in figure 1(d), and then show that the fourth velocity distribution 
can be found self-conslstently. By setting the reflected ton population 
to zero, and applying the Penrose (i960) criterion, they show that the 
plasma Is stable against small perturbations If the double-layer is long 
enough , and the transmitted Ion population is fast enough. 

In this paper, the double-layer Is analyzed using fluid theory. 

Section 2 treats the case of cold particles; thermal effects for the 
reflected particles are included in §j, and for the transmitted particles 
In §4. The emphasis is on finding a relationship between the magnitude of 
the potential step, and the parameters characterizing the adjoining plasmas. 

In § , qualitative statements are made about the length of the double- 
layer and a characteristic length useful as a quantitative measure is 
defined and discussed. 
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2 . COLD PLASMAS 

As a starting point, ail of the particles in figure 1(d) are taken to 
be at zero temperature. The "reflected" particles are actually stationary, 
and cannot penetrate the double-layer. Consequently, only transmitted 
particles need be considered In establishing the self-consistent potential. 
The density of the stationary particles will be determined so as to satisfy 
neutrality outside th*. double-layer. 

Electrons enter the double-layer from Region 1 with velocity u^ , 
density n^ and flux I~ e ( = n e l u e i^ • I° ns from Region 2 enter with 
velocity u t , density n^ and flux n i > u 1? ) * At a P olnt > x t 

within the double-layer, where the potential is <t> , the particle velo- 
cities and densities, determined from conservation of energy and particle 
f lux ,are 



where m^Cm^) is the electron (ion) mass, and e is the magnitude of the 
electronic charge. Substituting the densities frc i Eq . (l) in Poisson's 
equation yields. 



— (n -n ) 

C 0 e i 



Equation ( 2 ) can be used to define a function P(0) 


such that 


d g fl 

dx 2 


dP 
d© * 


Integrating, we obtain 


( 2 ) 


( 3 ) 
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where the constant of Integration Is chosen so that P(0) » 0 . Alterna- 
tively, P can be expressed as 


\ 


P(x) - .f 0 (x')E(x')dx' . (5) 

■'o 

The electric field Is assumed to vanish at the edges of the double-layer, 

E(o) =. e^) - 0 . (6) 

This ensures that the charge density Integrated across the double-layer 
vanishes. The function P Is thus constrained to satisfy 

P(x x ) - 0 , P(0 Q ) - 0 , (7) 


and Is negative within the double-layer. Since the Integrand In the last 
equality in Eq . ( ) Is the force between the double-layer and the particles, 
Eq . (7) can be viewed as a statement of mechanical equilibrium. 

Setting P(® ) = 0 In Eq . (4) yields 


2 mi r, 
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where M = m / m . Although Eq . (8) predicts a positive value for 0 

1 Q U 

within the ranges 
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the actual range is more restrictive. This can be seen by solving 

p(® 0 ) = 0 for r # /r t , 
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which varies monotonlcally from 
zero to Infinity. The rinses of 


1/2 

u ,/u to M as varies from 
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possible double-layer solutions are thus 



r 2 

r? > 

1 
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Langmuir's model (Langmuir 1929 i Block 1972) of the space-charge- 
11ml ted diode is contained In this formalism as a special case. Finite 
fluxes of particles enter the double-layer with zero velocity, so that 
the numerator of Eq . (8) vanishes. (The emissive powers of the cathode 
and the anode are assumed Infinite.) Therefore f 0 can be obtained 
only If the denominator also vanishes. This yields the Langmuir v 1929) 
condition 


r /r 


i 




1/2 


( 12 ) 


and leaves the potential difference Indeterminate. 

Equation (8) may be rewritten with the fluxes eliminated In favor 
of the particle densities, 


$ c 2 NO 


( 1 -NW ) ( N-l ) 

(i-Tor 


where we have Introduced the normalization 
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e* 
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e el 
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m u ' 
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e el 


( 13 ) 
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Equation (11) reduce* to the limits 

N > 1 > if Q , N < 1 < IT Q . (1^) 

These regions in (N,Q) space are shown shaded in figure 2. 

For neutrality in Plasmas 1 and 2, the reflected particle populations, 
normalized to n^ , must satisfy 



where n,, is the density of reflected Ions In Plasma 1 and n . Is 
ll e? 

the density of reflected electrons In Plasma 2. Since and 7]^ 

represent particle densities, they must not be negative. For N <. 1 , 

7| Is always positive, while for N > 1 , 7) e Is always positive. I’slng 
Eq . (13); and examining the non-trlvial cases yields 

N^Q > ^ (N < 1) ; JtQ < ~~ (N > 1) . (IT) 

The conditions of F.q. (17) are less stringent than those for existence 
of a double-layer, Kq. (15). It follows that for double-layer formation 
there must be reflected particles on both sides of the double-layer. 

Figure 3(a) shows the electrical potential, i(7.) , electric field, 

E (Z) , and charge density, p(7), for N ■ 1.1 and Q * 0.8. The normaliza- 
tion for the spatial coordinate, Z , is 



The electric field and charge density are normalized so that 

0 = dE/dZ = - d 2 $/d/. . 

Although the precise spatial variation of the electric potential 
can be determined only by Integration of Eq . (3); some symmetry pioperties 
can be deduced from the function II , 
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Tor the electric potential to be symmetric about 1 /2 , H must also be 
symmetric about { /? . This Is equivalent to requiring all odd derivatives 
to vanish at that point, 
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l-?m l-2m 

* (2m-i)!! 

l 

* * t-u* U* « 0 ) * 

V 2 

Q 
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(m odd). 


0 , 

( 20 ) 


This rearranges to 



/ o 

(2 


1) 


which Is satisfied by X 0 1 or jf 0 1 , where ! ** ® , only. 

The asymmetry of the electric potential can be gauged by determining the 

value of { at which the minimum of II occurs; this corresponds to 
M 

vanishing charge density and thus to the maximum electric field, 


jl n'q * K Q - + 1 

i -*n( nq-1) 


Hy expanding for weak (N ^ 1) and strong (ft O 
these cases are found to be nearly symmetrical, 


( 22 ) 


■v 1) double-layers, 


A 






6 

T— 

4 


D ___ 

(N - 1 + 6) , r - * (^Q . 1 + 6) . 


Contour plots of { in the (N,Q) plane are shown In figure *,and 

Indicate $ /$ , 1/2 over a broad range of N and Q . 

M U 

To summarize the foregoing results for cold particles, we note that 
the double-layer may be completely characterized by the densities and 
energies of the Ions and electrons streaming into It. The densities of 
the stationary, "reflected” particles, and the magnitude and spatial 
variation of the electric potential, are uniquely determined. 
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3. RI.'F LETTED PARTICLE TEMPERATURE NON-ZERO 

Thermal spread of the reflected particles will be described by assuming 
a Holt/.mann distribution within the double-luyer . The temperatures of the 
reflected ions and electrons, expressed In energy units, are T^ . and T^ 
respectively. Although this allows an exponentially small number of 
"reflected" particles to cross the double-layer, for T^ , , K< e<? 0 . 

this error will be negligible. 

The non-zero temperature of the reflected particles allows them to 
penetrate and partially neutralize the double-layer. The sharp boundaries 
associated with cold particles are thus smoothed out. The potential drop 
across the double-layer is also affected. 

The normalized Poisson equation is 



where ( T . . /m u' ) and 3 ( T _/m u are the normalized Ion and 

1\ He el I e \ e2 e ell 

electron temperatures. 

Equation (2U) is Integrated to find II , corrected for the non-zero 
reflected particle temperatures. Noting that Lq . (1- ) is still the proper 
ideu :if ication for t, and T , It is lound that 


n(*) - 1 - 


1/2 

( 1 + 2 *) 




JL 

( 1 ^ 2 * 0 ) 






( 25 ) 
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The self-consistent value, i , Is again determined from the require- 
ment 11(4, ) 0 . In tills case, 4 must be found numerically. For very 

low temperatures, an approximation to the difference, 6$ , between the 
true value, } , und that based on Bq . (13), 4^ , for cold particles, 1 b 

found to be 


64 *- o 


rr 

fj 

N(l+24 c ) 1/? - 1 

+ 3 

(>*S',r- 1 
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_*(l4S« c ) 1/S (l * | O' ^ -1. 
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7 77T 375 — 
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(26) 

At either edge of the rimible-1 nver, both T1 and d*Vd4 vnnlph; near the 
edge , II Is thus approximately 


d4 


4=4 


(*-0,* 0 ) • 


(27) 


p p 

Since 11 must be negative within the double-layer, d' II/d4" Is required 
to be negative at 4=0,4 From Eq . (2‘ ) the derivative Is 


d 2 n 


-,-(! + 24) 


-3/2 



16*1 Co-)) 


-3/2 




(28) 


"his reduces to 
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- 1/2 

N - (1 + 2 4 q ) 

h < 1 + 2 V" 3/2 


(*-» 0 ) • 
(29) 


The temperatures are thus bounded from above, JJ. < 1 and 3 • 0 . 

1 2 

In terms of unnormalized variables, these become T ^ < m e u e ^ and 
< m^u , which may be recognized as Bohm ( 19 ^ 9 ) conditions for 
collection of electrons and ions, respectively, through a sheath. 

The limits expressed by Eq . (29) must be compared with the physical 
requirement that the trapped particle densities not be negative. Consider- 
ing N as the dependent variable, and 4 as an independent variable, 

Eq . (29) may be rewritten as 
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(Q-3 e )(l + ?* 0 ) J 


These bounds on N are more stringent than those derived from the require- 
ment that the trapped particle densiti ; s not be negative. 


(l ♦ 2 * \ > » > 

\ Q 0 (1 + 
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As for cold particleu, a self-consistent solution for the double-layer 
requires trapped particles on both sides of the double-layer. 

Figure 3(b) shows 4(z) , E(z) and o(Z) for the values o' N and Q 
used in figure 3(®)> but with warm reflected particles. The penetration 
of the double-layer by the reflected particles produces the smooth variations 
of 4 , E and p . 

Examination of 11(4 ) for small 4 indicates the possible existence 
of a new root 11(4 ) = 0 , and thus a new double-layer solution (see figure 5 ), 
If the solution satisfies Eq . (jO), it is physically admissible. The catise 
of this modification may be identified by considering 11(4 ) as a function 
of 4 , and evaluating its derivative at 4=0 from Eq. (25). We obtain 


<*n(4 0 ) 

»o-° 


= 1 - N . 


This compares with dll(f )/d4 t = N-l t 4=0 when the reflected 

particles are cold, as can be seen from Eq . (19)* Each of the reflected 
species contributes a term 1-N , so that the slope of 11(4 ) is changed 
in sir.n by the inclusion of non-zero temperatures. 

If only one of the thermal terms is retained, the first derivative 
vanishes and the second derivative is found to be 


<rn(4 0 ) 


N-l N 1 


(% * °) , 


V 3 


(* t * °) 
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The value of 11(4 ) aa 4 approaches Infinity la unaffected by 
the finite temperatures: 

n (* 0 ) - ( 2» 0 ) l/2 '• (* 0 ~ *> • (34) 

For example, If the Initial alope of 11(f) Is negative [N> 1 , from 
Eq. ( 32 )] “nd the asymptotic value la positive [ IT Q > 1 , from Eq. ( 3 ^)]> 
there must be at least one root (in genera^ an odd number of roots) 11(4^) 3 0 . 
The analysis based on cold particles excluded any roots in this case. Similar 
reasoning shows that In the ranges limited by the cold particle analysis, 

A 

N £ 1 % ff Q , there will be either two roots or none (in general f an even 
number of roots). 

Figure > shows the behavior of $ as 3 p and 3^ are separately 

Increased from zero. The curves marked 3 and 3. are based on the 

e 1 

linear approximation, Eq. (2* ) . It Is seen that 4 decreases with increasing 

for N > 1 and with Increasing 3 for N < 1 . The upper bound on 
e 1 

the temperature for these cases comes from Eq . (29)# and not from the Bohm 
conditions, which are its high-potential limits. For certain values of 
N and Q , the second double-layer solution introduced by the non-zero 
temperature is physically admissible. As figure 7 shows, the two roots 
n($ Q ) = 0 coalesce at a finite 4^ and then disappear. 

Figure 8 shows (shaded) the regions In (N,Q) space within which double- 
layer solutions can be found for non-zero reflected particle temperatures. 

The boundary at N = 1(4 0) for cold particles is found to vary with 

temperature. The value of 4 along this boundary also varies, but is 
everywhere greater than zero. The >T Q = 1 boundary, with infinite 4 , 

is unaffected by the inclusion of temperature. 

To establish the validity of the double-layer solutions along the low 
potential boundary, an additional criterion should be investigated. If 

is not significantly greater than the reflected particle temperatures, 
a considerable fraction of the "reflected" particles will be transmitted 
through the double-layer, thus invalidating the model. In fact, all points 
along the boundary are found to satisfy 4 ^ > 3 . For reference, the line 
4 = 3 is shown in figure 8. 
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4 . TRANSMITTED PARTICLE TEMPERATURE NON-ZERO 


If the transmitted particle velocity distributions have narrow widths 

about their streaming velocities, a fluid model can be used to describe 

the particle dynamics. Denoting; the electron temperature by T , with 

e 

value T^ within Plasma 1, and the Ion temperature by T , with value 
T^ within Plasma 2 , the momentum transfer equations are 



Assuming that the particles are accelerated adlabatically, 


dT 

T 


(V-l) 


dn 

n 



( 36 ) 


where V 3 for one-dimensional adiabatic acceleration, the thermal terms 
on the right-hand side of Eq . ( 3 ‘ ) are 


1 d_ 
n dx 


(nT) 


3 dT 
2 dx 



( 37 ) 


where T and n are evaluated at some reference position. Using the 
continuity of particle flux for each species, 


u 


un 

"n > 


( 38 ) 


und Eq . ( 37 ), Eq . ( ) can be written as a conservation principle, 
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( 39 ) 
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With the normalization given in $2, the conserved quantities become 


\ O V* Q 3 

— - * ♦ “ T _>T . const, + . const, (40) 

2tf 2 e e 2Nj 2 W 


where the normalized temperatures and densities are 
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(41) 


The constants in Eq . (40) are found by evaluating * e left-hand sides in 
Plasmas 1 and 2 for the electrons and ions, respectively. The expressions 
in (40) can be solved to yield 



24 + 3% + 1 - [(24 + jr e + l) 2 

6 t 

e 



(42) 


1^(4) = N 2 


2(4 ( -4) + 3 T t + Q - 

(2(4 ( -4) + 3 T t + Q)‘ - 12 T^j 

1/2' 

_ 


. 


\ 


These expressions are valid only for 

3 T e < 1 , 3 T 1 < <3 • (43) 


For higher temperatures, the fluid approach breaks down, as discussed in 
Appendix A. 

The densities given by Eq . (42) can be used in Poisson's equation to 
generate a new function n . The requirement of neutrality outside the 
double-layer, i.e. that dll/d? vanish at 4=0,4, yields appropriate 
values for the reflected particle densities, 


71. = 1 - N 


24 ( + 3T t + O - [(24 0 + 3T t + o)‘ - l2T t Q 


1/2 


1/2 


6 T . 


n = N 


2 V 1 - K* K* 1)g - **%] 


1/2 


1/2 


(44) 
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The function II can be divided into four component*, one for each of 
the transmitted and reflected classes of electrons and ions, 



The self-consistent value of the potential step, 4 , must again 

.'e found numerically. Noting that dC(4 )/d4 1 - N at 4 = 0 , and 

1/2 1/2 0 0 0 
that H -• (24 ) ' ( NQ - 1) as 4 m > we 8ee that the gross behavior 

of fl(. ) is unaffected by the inclusion of the transmitted particle 

temperatures. No new double-layer solutions are to be expected, though 

their location and physical admissibility may be modified. 

The change, 64 , between the vtlue predicted for cold particles, 4^ , 

and the true 4, can he approximated for small t and t, hv 
0 e i 7 



Equations (°^) and (h£) are additive when all four species have finite 
temperatures . 

Figure c shows the variation of 4 as and T ^ are separately 

increased from zero. The curves marked and T are based or the 
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linear approximation, Ec . (1*0). The qualitative effects are similar to 

those for non-zero 3 and 3 though 4 Is not double-valued and has 

e 1 

no low potential cutoff. The condition d^Il/dt < 0 for 4 ■ 0 , 4 ( 

becomes 


dN 
c 

d4 


\ l V 4 \ \ /»\ "i 

- % exp (- -vr 5: exp U) 


(♦ > 0,0 , (UT) 


where dN '6$ and dNVdJ are evaluated 1 rom 1>< . fU?), 
e 1 


dN 

e 

d$ 


^1 

d t 


24 + 3t 

e 

+ 1 - 

(24 ♦ *■ I)' 

- I£f 

0 m 

1/2 

6t 

e 

’(24 + 

3t * i) : - lr t 

e e 



i -;) • ■f 1 ■ o - ( : ^ *q) - U 

6T i[( 2( V > ) " 3T 1 * °) 2 - ^i 9 ] 


(U8) 


-.1/2 


Since dN /d4 Is negative and dN /d4 Is positive, l.e. electrons 
e 1 

become less dense and ions more dense In moving to higher potential, T) 

'e 

and T)j have positive lower bounds. In the limit of strong double-layers, 

4 — t> , Eq. (47) reduces to 


+ 3T e < 1 (4=0) , 3 e + 3T i < 9 (*"V • 

or, In terms of unnormalized variables 


T il + 


el 


m u' , 
e el 


e2 


3T 12 < m i 12 


(*9) 

(50) 


These upper bounds apply only In the Infinite potential limit; for finite 
potentials, Eq . (U7) is more restrictive and is thus the appropriate 
criterion to be applied. 

If the reflected particles are cold, dll/d4 is discontinuous at 
4 0,4 even though the transmitted particles are warm. Equation (1*7) 

is then replaced by the requirement that the density of the reflected 
particles, [*f] , 71 e from Eq . (1*1*)], not be negative. 


; 
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. DOUBLD- LAYER DIMENSIONS 


Having determined how the potential drop, I , of the double-layer 

V 

oependa on the physical parameters of the adjoining plasmas, the remaining 
question accessible to this analysis is how the spatial extent of the double- 
layer is determined. 

In determining the distance required for the potential to vary from 
0 to , two distinct cases can be distinguished by making a power- 
series expansion for T. . If the reflected particles are cold, the charge 
density, p( a dH/di) , will be discontinuous at the double-layer edge. 

We thus have 

n « - a* , 



- 

d i 


as a , 


(t > 0 ) 


( 51 ) 


which can be solved to yield 

* s jr a z 2 (e > o) . 


( 52 ) 


From Eq . (‘>2), we see that a finite distance is required for the potential 
to vary from zero to a given non-zero value. 

If the reflected particles are warm there will be no discontinuity 
In p at the double-layer edge. We thus have 


II as - , 

(* > 0 ) ( 53 ) 



which yields 

{ s\ exp|(2g)^ zj (* l 0) . 


( 5«0 


The exponential dependence in Eq . (54) indicates that the potential approaches 
zero for any finite Z , but attains it only asymptotically at Z * - * . 

The behavior at i - 0 is determined by the reflected ions in 
Plasma 1; similar arguments apply at $ = 4 , with the behavior 

dependent on the reflected electrons in Plasma 2. Therefore, with warm 
reflected part icl es, the double-layer is infinitely long, i.e., electrical 
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neutrality obtains only at 7. - ± » , even though most of the potential 
variation occurs over a finite length. A general method of defining the 
length of a double-layer Is thus required. 

Andrews c* Allen (l9fl) dellne the double-layer length as the distance 
over which the potential varies from a tenth of the trapped Ion temperature 
on the low potential side to within a tenth of the trapped electron tempera- 
ture on the high potential side of Its asymptotic value. Hasan et ter Haar 
(I 978 ) measure between points on both sides of the double-layer where the 
potential Is within a tenth of the trapped electron temperature of Its 
asymptotic value. Knorr Goertz (197*0 assume a hyperbolic tangent form 
for their analysis of a double-layer, 4 = $ r tanh(x/f) , and therefore have 
the built-in length scale, ? . 

We choose to define the characteristic length, L , of a double-layer 
as the length over which the constant electric field evaluated at 4 q/ 2 
would have to exist in order to produce the potential difference 4 , 

as Indicated by the dashed lines In figure 3> 



While this definition is no less arbitrary than previous definitions, it 
commends itself as requiring less numerical integration, and defines the 
region of strong electric field, over which most of the potential varia- 
tion occurs. 

Applying Eq . (‘ ; ) to the case of cold particles, the length of the 
double-layer is found to be 


L = 


(- 2 [l - (l . » Q ) 1/2 * NQ(1 * §■ »J 1/2 - N9(l - ^-) 1/2 ]) 1/S 


(56) 


When the double-layer is very strong, 4 Q - » [ N^Q 1 , see Eq . (13)], this 
reduces to 


L =- 0.92 4 


3A 


(57) 


or, in unnormalized variables, 
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i 



0 3A 

o 

T7F 


( en ol U el ) 


This differs by a factor of order unity from Child's Law for a diode, 
where en^u^ , 18 the electron current density. This is consistent with 

the assumption of a large potential difference, i.e. that the partlcleB 
enter the double-layer with negligible kinetic energy relative to that 
with which they leave, ke note further thut l^Q - 1 is the normalized 
form of the Langmuir condition, Eq . (12). 

For weak layers, 4 -* 0(N -• 1) , Eq. (56) reduces to 

1 “ 2 (&) • (59) 

so that the length of very weak double-layers varies smoothly from L a 0 to 
L 3 2 as Q varies from zero to infinity. 

Figure 9 shows the dependence of the length of the double-layer on 
the parameters describing the adjoining plasmas. In those cases where an 
increase in one of the temperatures causes a decrease in the potential 
difference across the double-layer, there is an associated decrease in its 
length as well. 

It is important to note that distances are not normalized to the 
electron or ion Debye length, but to the distance a transmitted electron 
travels in a plasma period, before it enters the layer. The scaling is 
thus independent of temperature. It follows that no difficulty arises in 
treating the case of cold particles, where the double-layer extends over 
an infinite number of Debye lengths. 


\ 
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t DISCUSSION 

it Is Instructive to apply our model to conditions characteristic of 
laboratory and space plasmas that may sustain double-layers. Although 
the physical parameters, l.e. density, temperature, and drift velocity, 
differ by orders of magnitude, the normalized variables may be similar. 

For the doub le-luyers observed by Quon a* Wong ( 1976) .we make the 
rough estimates 

N 1.2 , Q a; 0.5 , T 0.2 , T. * 0.1 , 3 * 0.1 , 3. * 0.3 , (60) 

v 1 6 1 

which produce a double-layer of magnitude and length 4 * 4.9 , L ■ 5.6 

,8-3 0 

For a density of 10 cm , and electron energy of 1 eV in Plasma 1, this 

corresponds to a potential difference of 20 V with a characteristic distance 
of 0.3 cm. This Is somewhat stronger than the 3-1 V and sharper than 
the 3- cm reported, but of the right order of magnitude. 

As we have shown, comparison of the admissibility criteria with the 
requirement that the reflected particle densities not be negative indi- 
cates that reflected particles must be present on both sides of the double- 
layer. In their experiment, Quon a, Wong were able to remove the reflected 
Ion population by biasing a grid appropriately. They report that, without 
the reflected Ions, they could not produce a double-layer. 

Although we lack a complete set of measurements made during an active 

2 

aurora, we assume a current of 1 uA/m may be carried by a 100 eV electron 
beam above the double-layer. For parameters 

N 1.* , Q = 0.4, t = 0.1 , t 3 0.1 » ** e * 0.1 » 3 * 0.1* 

( 61 ) 

(compare with values used by Swift 1976), a double-layer of magnitude and 
length, 4 = 10.3 , L » 7 .4 , results. This corresponds to a 2 kV 

potential in a 0.7 km long step, which requires an average electric field 
of 2.8 V/m. This field is a factor of five greater than that measured by 
satellite (Mozer et al. 1977 ), but the measurements may not have been made 
In the center of the double-layer. If the current above the double-layer 
is carried by 1 keV electrons, the strength of the double-layer increases 
to 20 kV, while the electric field increases to only .0 V/m. These 
potentials are of the right magnitude to produce the observations of 
high-energy electron precipitation (see, e.g. Shawhan et^ a_l. 197 > a °d 
references therein) . 




For self -con* latency the relationship* of Lqs. (29) and (1*7) between 
the particle energies and temperatures have been developed. If tlieae 
condition*) are not satisfied outside the double-layer, u "pre-sheath" 

(Schott 1 * ) may be postulated thut would accelerate *.ne particles suffi- 

ciently. The excitation of pulses at the Ion acoustic speed, as seen In 
some positive column discharge (liable at Torven 197*0 ma V originate In 
such a pre-sheath. 

The reflected particles with non-zero temperature have been assumed 
to obey Boltzmann's law within the double-layer. This can result either 
from an Isothermal equation of state In a fluid theory approach .or from 
a Maxwellian velocity distribution In a kinetic theory approach. Non-zero 
temperature of the transmitted particles has been handled in a fluid theory 
approach, with an adiabatic equation of state. However, as Is shown In 
Appendix U, the same results can be obtained from a kinetic theory approach 
with vaterhag distrilotion functions; results for other distribution 
functions are considered, and compared with the results of fluid theory. 

One Important aspect of the double-laver that Is not treated within 
this time-independent analysis Is that of stability. The Penrose (1960) 
criterior for determining stalllltv applies only to homogeneous systems. 

It may thus be possible to find velocity distribution functions that are 
Penrose stable, but which are unstable due to the inhomogeneit v (Wahlberg 
1977). .'Iso stability, as determined from the Penrose criterion, depends 
critically on the distribution functions. As a solution of Poisson's 
equation, the double-layer is dependent only on t lie density of the particles, 
l.e. the integral of the distribution functions. Furthermore, there Is 
evidence from both laboratory (Hal if & Torvort 1974; Armstrong 1975; Quon & 
Wong 1976) and romnuter exnerinents (Hubbard 6 Joyce 1°78) that double- 
layers can co-exist with some (ostab 1 1 ity . 

The analysis presented here stun es the self-consistency of the double- 
layer as an independent entity. However, since double-layers form a part 
of a larger circuit, whether in the laboratory or in space, there are 
further problems relating to consistency with the external circuit that 
will reouire investigation in future work. 
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APPENDIX A: LIMITS OP FLUID THFOPY 


For y i 1 In (3*), (37) Is rewritten 


1 iL 

n dx 


(nT) 


„ _L_ 41 . _1_ d_ 

A 

T 

[Cl 

Y-l dx y-1 dx 


I'J 


(A.I) 


The Isothermal equation of state, y « 1 , will he treated separately, 
below. 

The conserved quantity analogous to (40) Is 


2n* 




i 

■ h 




(A. 2) 


where q Is the charge of the particle, u Is the drift velocity, T 
Is the temperature, and f Is the potential at the reference position. 
Rearranging (A. 2) yields 


F 

y 


* 


n) 

*? 

A 

n 

T 

■ u' 

— 

n- 


n 



Y-l i" |n. 


u‘ ♦ ^ (♦ •#) 



t 

m 


(A. 3) 


The function F (n/n) is niotted In figure A. 1(a), and shows a minimum at 


f nl 

f 


M “ 

'n- 1 m 

[vT] 


(A. 4) 


s 


As qf decreases, the density must Increase or decrease as Indicated in 
figure A.I in order to satisfy (A. 3). However, decreasing a: represents 
acceleration which, to conserve particle number, must he accompanied by 
a decrease in the density, n . Values of n/n greater than (n/n) 

m 

are therefore nonphysical. Thus, evaluating (A. 4) at the reference position 
yields 

y T < mu - . (A. 5) 

Applying (A. 5) to the electrons entering the double-layer from their 
reference point in Plasma 1, or to the ions entering from their reference 
point in Plasma 2, and normalizing as in (41), yields 
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This is the result of (\.4) for >“1 , and therefore implies the temperature 
linitations of (A.b) afiain. 
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APPENDIX Bs COMPARISON OF FLUID AND KIVETIC MODELS 
FOR TRANSMITTED PARTICLES 

Instead of using a fluid model for the transmitted particles, and 
then assuming an equation of state, it is possible to assume a form for 
the velocity distribution function, f(u) , before the particles enter the 
double-layer, and determine its dependence on the potential, $ , by the 
Vlasov eauatlon. Taking f(u) as time-independent, the Vlasov equation, 
in one dimension, is satisfied for any distribution that is a function of 
the total energy; 


f(u,<t) 


- f 

+ 

u 2 + 2oi 



m 



(B. 1 ) 


where m and q are the mass and charge of the particle. The spatial 
dependence of f enters through the potential, $ . The sign of the argu- 
ment in (B.l) is the sign of tbt velocity. 

A simple choice for the distribution function at <t m 0 is a waterbag 
with all pa. tides forward-going, [see figure B.l (a)] 


"o 

l 6 U t 

, . ,1/2 
|U-u D |< 3 u t . 


(° 

1 1 > 1 / 2 
|u-u D | >3 u t . 

(B. 2) 


The Identifications of n^ as the density, u^ as the drift velocity, 
and u ^ as the thermal velocity, are made by calculating the appropriate 
moments of the distribution. 

Normalizing the density and potent ial r and defining a temperature con- 
sistent with §4, 


4> * , T 


mu: 


(B.3) 


the density as a function of 4> is 

.1/2 (( r , 2 


N($) 


6t 


1 72 


1 + (3t) 1 /2 | + 

J 


2 $ 


1/2 


1 - (3t ) 


1/2 


1/2 


+ 2 < 


(B.4) 
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vhich is equivalent to (42) for electrons. The assumption of all forward- 
going particles is 


3i < 1 , 


<*.5) 


which is (43) for the use of fluid theory. Thus as far as the analysis 
in §4 is concerned, the waterbag and fluid models are identical. 

If the velocity distribution is allowed to have some backward-going 
particles, as in figure B.l(b), the density, normalized according to (B.3) 
is 


N(0 





" 1^/2 i >2 

+ 2t\ -(8t) lU 


(B.6) 

Although the mathematics does not allow the general case of a drifting 
Maxwellian to be treated analytically, two particular cases can be studied. 
If the distribution is as in figure B.l(c), a Maxwellian with a large 
drif t-to-thermal velocity ratio, and thus with a negligible number of 
backward-going particles, the distribution function is 


f(u.O) 


u t ( 2v ) 


1/2 


exp 


2u 2 


(B. 7) 


To second order in t , the normalized density is 


»(*) 


1 

(1 + 2 *) 


T72 


1 


3$t 

(1 + 2 $) ? 



(3* -2) St 2 
(1 + 2i) 



(B.8) 


This agrees to first order with the small-T expansion of N for the 
f luid /waterhag model, the second order term of (B.8) being larger by a 
factor of 5/3 than in the f luid/waterbag model. 

For the forward-going half of a Maxwellian at <f“0, illustrated in 
figure B.l(d), the distribution function, drift velocity and thermal 
velocity are 


26 


<u > 0) 



(u < 0) , 


U D " U 0h] 


1/2 


= 0.80 u 


0 ’ 


f tt- 2'I 1/2 * 0.60 u. . 
u ■ u_ 0 

t 0 { v J 


(B.9) 


It can be seen from (B.9) that the normalized temperature is a constant, 
•t . ( ft — 2 ) / 2 = 0.57. The normalized density, as the particles are accele- 
rated, is 


n( 4 ) 


exp 



1 



(B. 10) 


Figure B.2 shows the dependence of the density on the potential 
for the various distribution functions considered here. Foi all but 
r ■ 0.3, the f luid/waterbag and laxweilian curves totally overlap; the 
Maxwellian is the lower of the t = 0.3 curves. Since Poisson's eauation 
depends only on the density, and not on the detailed shape of the distri- 
bution, the fluid /waterbag model may be used for a Maxwellian with a large 
dri f t-to-therral velocity ratio, but not for a half-Maxwellian. 
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FIG. 1. Various double-layer models. 

(a) Transition region between two plasmas with a potential 
difference. 

(b) Space-charge-limited diode. (Particles enter with zero 
velocity, but with non-zero flux.) 

(c) Space-charge-limited cathode (Cathodic electrons enter with 
zero velocity, but with non-zero flux. Plasma electrons 
are reflected by the double-layer; ions are transmitted 
through it.) 

(d) Double-layer between two plasmas. (Four species of particles 
are reflected and transmitted.) 

(e) Double-layer between two plasmas. (Ions from Plasma 1 are 

transmitted by the double-layer. Only one species of particles 
is reflected: electrons in Plasma 2.) 
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FIG. 5. Double-layer s lutlons for non-zero reflected particle 
temperatures. 

(a) Modification of cold plasma solution. 

(b) Additional solution. 
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FIG. 9. Double-layer length. 

(a) Cold plasma theory. 
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FIG. 9. 
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(Contd.) Double-layer length, 
(c) Macroscopic plasma theory. 
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FIG. A.l. Variation of functions F and F 

y 1 

(Arrows correspond to q0 decreasing) . 

(a) Adiabatic. 

(b) Isothermal. 
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FIG. B.l. Velocity distribution functions at q<f> ■ 0 and q$ < 0. 

(a) Waterbag with all particles forward-going. 

(b) Waterbag with forward- and backward-going particles. 

(c) Maxwellian with large drif t-to-thermal velocity ratio. 

(d) Half-Maxwellian . 








FIG. B.2. Variation of N(<t>) for the velocity distribution 
functions of figure B.l. 


